Structural relaxation in a system of dumbbell molecules 
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The interaction-site-density-fluctuation correlators, the dipole-relaxation functions, and the mean- 
squared displacements of a system of symmetric dumbbells of fused hard spheres are calculated for 
two representative elongations of the molecules within the mode-coupling theory for the evolution of 
glassy dynamics. For large elongations, universal relaxation laws for states near the glass transition 
are valid for parameters and time intervals similar to the ones found for the hard-sphere system. 
Rotation-translation coupling leads to an enlarged crossover interval for the mean-squared displace- 
ment of the constituent atoms between the end of the von Schweidler regime and the beginning 
of the diffusion process. For small elongations, the superposition principle for the reorientational 
Q-process is violated for parameters and time intervals of interest for data analysis, and there is 
a strong breaking of the coupling of the a-relaxation scale for the diffusion process with that for 
representative density fluctuations and for dipole reorientations. 

PACS numbers: 64.70.Pf, 61.20.Lc, 61.25.Em 



I. INTRODUCTION 

Recently, the mode-coupling theory (MCT) for the 
evolution of glassy dynamics in systems of spherical 
particles has been extended to a theory for systems of 
molecules. The fluctuations of the interaction-site densi- 
ties have been used as the basic variables to describe the 
structure of the system. As a result, the known scalar 
MCT equations for the density fluctuations in simple sys- 
tems have been generalized to n-by-n matrix equations 
for the interaction-site-density correlators, where n de- 
notes the number of atoms forming the molecule. The 
theory was applied to calculate the liquid-glass phase di- 
agram and to evaluate the glass form factors for a hard- 
dumbbell system (HDS) |J. In the following, the pre- 
ceding work shall be continued by evaluating the con- 
ventional time-dependent correlation functions near the 
liquid-glass transition. It is the goal to examine the range 
of validity of the universal relaxation laws and to identify 
features of the glassy dynamics which are characteristic 
for molecular as opposed to atomic systems. 

The dumbbells to be studied consist of two equal hard 
spheres of diameter d which are fused so that there is 
a distance (d, < ( < 1, between their centers. The 
system's equilibrium structure for density p is specified 
by two control parameters: the elongation parameter £ 
and the packing fraction tp = ^pd 3 (l + §C ~ |C 3 )- The 
liquid-glass transition points f c (C) form a non-monotonic 
<^ c (£)-versus-£ curve in the (^-^-control-parameter plane 
with a maximum for £ near 0.43. There is a second transi- 
tion line <£a(0 within the glass phase tp > tp c (()- It sepa- 
rates a plastic-glass phase for <p c {0 < <P < Va(£)i where 
dipole motion is ergodic, from a glass for tp > <pa(C,)i 
where also the molecular axes are arrested in a dis- 
ordered array. The second line <£>a(C) terminates at 
C = C = 0.345, where <p A ((c) = ^c(Cc) {of. Fig. 1 of 
Rcf. [y). In the present paper, it shall be shown that 
there are two scenarios for the liquid-glass transition dy- 
namics for £ > Q c . The first one, to be demonstrated for 



C = 1.0, deals with strong steric hindrance for reorien- 
tational motion. For this case, all universal relaxation 
laws hold within similar parameters and time intervals 
as found for the hard-sphere system (HSS) |§, ||- As 
a new feature, there appears a very large crossover in- 
terval for the a-process of the constituent atom's mean- 
squared displacement for times between the end of von 
Schweidler's power law and the beginning of the diffu- 
sion regime. The other scenario, to be shown for £ = 0.4, 
deals with weak steric hindrance for reorientational mo- 
tion. The universal laws for the reorientational a-process 
are restricted to such narrow tp — tp c (Q intervals that 
they are practically irrelevant for the interpretation of 
data obtained by molecular-dynamics simulations or by 
presently used spectrometers. 

The paper is organized as follows. In Sec. II A, the 
MCT equations of motion for the coherent and inco- 
herent density corre lation functions for the HDS are 
listed. Section II B contains known universal laws for 



the MCT glass transition. The next section presents the 
new resul ts for representative density-fluctuati on cor rela- 
tors (Sec. [II A ), for the dipole dynam ics (S ec. Ill B ), for 
the mean-squared displa ceme nts (Sec. Ill C| ), and for the 
a-relaxation scales (Sec. [II D). The findings are summa- 
rized in Sec. IV. 



II. MCT EQUATIONS FOR THE HDS 
A. Equations of motion for the density correlators 

In this subsection, the basic equations for the system 
of symmetric hard dumbbells are noted. They have been 
derived in Ref. [j], and their solutions underlie all results 
to be discussed in the present paper. 

If r", a = A or B, denote the interaction-site cen- 
ters of the ith AB dumbbell molecule, the interaction- 
site-density fluctuations for wave vector q read pi — 
53i ex P(*9 ' ?t)- Similarly, the tagged-molecule-density 
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fluctuations read p~ s = exp(iq ■ r s a ) with r" denot- 
ing the interaction-site positions of the tagged molecule. 
It is convenient to transform to total number den- 
sities p~ = (p~ + p?)/^/2 and "charge" densities 

p? = (p4 — p§)/^/2. Similar definitions are used for 
the tagged-molecule densities. The top-down symme- 
try of the molecules implies vanishing cross correlations 
0#(*)*pf (0)) = {p^MpIM) = and also the iden- 
tity (pf(*)*pf(o))/(|pf I 2 ) = (pl s (t)*PlM)K\pl s \ 2 )- 

There is only one independent coherent density correla- 
tor, which shall be used as normalized function 

4$V) = (f$V)*f%(P))/(\P3 1 2 )- (!) 

There are two incoherent density correlators 



q.s 



x = N, Z. (2) 



Because of rotational invariance, the density correlators 
depend on the wave- vector modulus q = | q | only. The 
short-time asymptotics of the correlators read: <j)q (t) = 
1 - ±(^i) 2 + 0(|<| 3 ) and 4> x q , s (t) = 1 - ^jf + 
0(\t\ 3 ). The characteristic frequencies are given by the 
formulas: (f!^) 2 = g 2 { V 2 [1 + j (<?C<*)] + (v 2 R C 2 d 2 /6)[l - 
Jo(qCd)-j2(qCd)]}/S q v and (fij J 2 = q 2 {v 2 T [l±j (q(d)} + 
(v%( 2 d 2 /6){l T 3o(q(d) T 32(q(d)]}/w* q for x = N,Z. 
Here ji denotes the spherical Bessel function for index 
£, abbreviates the total-density structure factor, and 
Wg = 1 ± jo(q(d) for x = N,Z are the intramolecular 

structure factors. Furthermore, vt = yJksT /2m denotes 
the thermal velocity for translation of the molecule of 
atomic masses m at temperature T, and vr — (2/C,d)vT 
is the thermal velocity for rotations. The Zwanzig-Mori 
equations of motion ||] for the specified correlation func- 
tions are: 

d 2 ^(t) + (^f^(t) 



+ (^) 2 f dt'm^{t~t')d t ,^{t') = 1 (3) 
Jo 



+ (^, s ) 2 / dt'ml s (t-t')d t ,<j>l s (t') = 0, (4) 
Jo 

for x = N, Z. All complications of the dynamics of the 
many-particle problem are hidden in the relaxation ker- 
nels rriq (t) and m x q s (t). 

Within MCT, the relaxation kernel (t) is expressed 
as functional of the density correlators: 



m"(t)=jW(i)]. 



(5a) 



The mode-coupling functional reads J 7 ** [/] = 

| J dkV N (q;k,p)fkf p where k + p = q, and the 
positive coupling vertices are given by the density p, the 
structure factor S q , and the direct correlation function 
Cq (cf. Eq. (20b) of Ref. Wave- vector integrals are 



converted into discrete sums by introducing some upper 
cutoff q* and using a grid of M equally spaced values for 
the modulus: qd = h/2,3h/2, ■■■ , (q*d - h/2). Thus, q 
can be considered as a label for an array of M values. 
Equation (||) then represents a set of M equations which 
are coupled by the polynomial 



k,p 



Vq,kp fk fp- 



(5b) 



It is an elementary task to express the M 3 positive coef- 
ficients Vq.k-p in terms of the V N (q;k,p) (cf. Eq. (7) of 
Ref. ||). Similarly, one derives functionals for the kernels 
for the tagged-molecule functions: 



x = N,Z. 



(6a) 



Again, the functionals Tq s are given in terms of the equi- 
librium structure functions (cf. Eq. (23) of Ref. Q). After 
the discretization, the mode-coupling polynomial has the 
form 



^g,fep 

kyP 



fk,s fp> 



N, Z. 



(6b) 



The following work is done for a cutoff q*d = 40 and 
AI = 100 wave-vector moduli. The structure factors 
are evaluated within the reference-interaction-site-model 
(RISM) theory @, § §. In Ref. 0, the details of the static 
correlation functions have been discussed. Equations (||) 
and (^|) are closed non-linear integro-differential equa- 
tions to calculate the M correlators <j)g (t). Using these 
correlators as input, Eqs. (||) and @ are closed equations 
to evaluate the M correlators <f>g ta (t) for x = N and Z. 
The mathematical structure of Eqs. (^) and (|J) is identi- 
cal to that discussed earlier for the density correlators of 
the HSS Q]. The differences are merely the values of the 
frequencies fl^ and the values for the coupling constants 
Vq.kp- A similar statement holds for Eqs. (Q) and (|^), 
which are the analogues of the equations for the motion 
of a sphere in the HSS [3|. Therefore, all general results 
discussed previously [§, 3| hold for the present theory as 
well. 

In the rest of the paper, the diameter of the constituent 
atoms shall be chosen as unit of length, d = 1 , and the 
unit of time is chosen so, that Vt = 1. 



B. Universal laws 

This subsection compiles the universal laws for the dy- 
namics near the liquid-glass transition. They will be used 
in the following Sec. Ill to analyze the numerical solutions 



of the equations of motion. The derivation of these laws 
is discussed comprehensively in Refs. || and 0, where also 
the earlier literature on this subject is cited. 

From the mode-coupling functional in Eq. (|5t|), 
one calculates an Af-by-M matrix C q k = 
{dF?[f N ]/df»Kl - f k N ) 2 . Here jf = <j$(t - oc) 
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denote the nonergodicity parameters for the density 
fluctuations of the glass states. This matrix has a 
non-degenerate maximum eigenvalue E < 1. The 
transition is characterized by E c = 1. Here and in the 
following, the superscript c indicates that the quantity 
is evaluated for ip = ip c {0- Let e and e denote the right 
and left eigenvectors, respectively, at the critical points: 
Y^kCqk e k = e g , J2 q e q Cq k = e k . The eigenvectors are 
fixed uniquely by requiring e q > 0, e q > 0, J2 q &q e q = 1; 
and ^2 q e q (l — f q c )e q e q = 1. These eigenvectors are 
obtained as a byproduct of the numerical determination 
of ifdO described in Ref. |l|. They are used to evaluate 
the critical amplitude 

< = (l-0 2 e 9 , (7) 
and the exponent parameter A, 1/2 < A < 1, 



A = 5 E «« {^"[/^ / a O <• ( 8 ) 

q,fc,p 

Furthermore, a smooth function <r of the control param- 
eters C and ip is defined by 

<* = E ^ Wt/^i - ^fi/^]}- ( 9a ) 



In a leading expansion for small tp — </? c (C)i one can write 



° = Ce, e=(^-^(C))M(C)- 



(9b) 



Glass states are characterized by a > 0, liquid states by 
er < 0, and a — specifies the transition. The nonergod- 
icity parameter exhibits a square-root singularity. In a 
leading-order expansion for small e, one gets 



f? = f? D + 



Nc 
1 



A), a > 0, 



0. (10) 



At the transition, the correlator exhibits a power-law de- 
cay which is specified by the critical exponent a, < a < 
1/2. In a leading-order expansion in (l/i) a , one gets 



4>"{t) = f q 



h™{t /t) a , (7 = 0, (t/to)->00. (11) 



Here, to is a time scale for the transient dynamics. 
The exponent a is determined from the equation r(l — 
a) 2 /r(l — 2a) = A, where T denotes the gamma function. 

Let us consider the correlator <pY(t) = 
{Y(t)*Y(0))/(\Y\ 2 ) of some variable Y coupling to 
the density fluctuations. Its nonergodicity parameter 
fy = 4>v{t ~ * oo ) obeys an equation analogous to 
Eq. ©: f Y = f Y + hy a J (1 A) + 0(a). The 
critical nonergodicity parameter f Y > and the critical 
amplitude hy > are equilibrium quantities to be cal- 
culated from the relevant mode-coupling functionals at 
the critical point ip = <p c (0- ^ ^ = P x qs -> ^ nc correlator 
4>y{t) refers to the tagged-molecule-density fluctuations, 
4> qs (t). Their nonergodicity parameters fy = f qs have 



been discussed in Ref. |l|, and the explicit formulas for 
the evaluation of hy = h q s can be inferred from Ref. ||. 
For small values of e, there is a large time interval, 
where 4>yit) is close to f Y . Solving the equations of 
motion asymptotically for this plateau regime, one gets 
in leading order in the small quantities 4>y(t) — f Y the 
factorization theorem: 



<f>Y(t)-fy = hyG(t). 



(12) 



The function G(t) is the same for all variables Y. It de- 
scribes the complete dependence on time and on control 
parameters via the first scaling law: 



G{t) 



\<j\g±(t/t a ), (7^0. 



Here 



t0 I \<T | 



l/2«, 



(13) 



(14) 



is the first critical time scale. The master functions g±(t) 
are determined by A; they can easily be evaluated numer- 
ically @. One gets g±(i -> 0) = l/t a , so that Eq. (0) 
for Y — p^i is reproduced for fixed large t and a tending 

to zero. Since g + (i — > oo) = l/Vl — A, also Eq. ( [To| ) is 
reproduced. 

One finds for the large rescaled time t: <7_(t — > oo) = 
-Bt b + 0(l/i b ). The anomalous exponent b, < b < 1, 
which is called the von Schweidler exponent, is to be 
calculated from A via the equation T(l + 6) 2 /r(l + 26) = 
A. The constant B is of order unity, and is also fixed by 
A 0. Substituting this result into Eqs. jl2] ) and (|T3|), one 
obtains von Schweidler's law for the decay of the liquid 
correlator below the plateau f Y : 



-{t) = fy-hy{t/t' a )\ * CT «t, 



-Q. (15) 



The control-parameter dependence is given via the sec- 
ond critical time scale t' : 



t'^toB- 1 ^ /|ap, 7 = (l/2o) + (1/26). 



(16) 



The leading corrections to the preceding formulas for 
fixed t — t/t a are proportional to \a\. They are speci- 
fied by two correction amplitudes and additional scaling 
functions h±(t/t a ) [^J. The dynamical process described 
by the cited results is called /3-process. The /3-correlator 
G(t) describes in leading order the decay of the corre- 
lator towards the plateau value f Y within the interval 
to ^ t <^ t a . The glass correlators arrest at fy for 
t 3> t a - In the limit a — > —0, all correlators cross their 
plateau f Y at the same time r^, given by 



Tfj = t- t c 



(17) 



Here, the number £_ is fixed by A via g^(tJ) = 0. 

The decay of <fiy(t) below the plateau f Y is called the 
a-process for variable Y . For this process, there holds 
the second scaling law in leading order for a — > — 0: 



M*) = M*)> t = t/t' a , t CT <t, 



(18) 
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which is also referred to as the superposition principle. 
The control-parameter-independent shape function cf>y (t) 
is to be evaluated from the mode-coupling functionals at 
the critical point. For short rescaled times t, one gets 



-(*) = n 



hyt" 



0(t 2b ), so that Eq. (jT|j is repro- 



duced. The ranges of applicability of the first and the 
second scaling laws overlap; both scaling laws yield von 
Schweidler's law for t a <C t <C t' a . 

The superposition principle implies coupling of the a- 
relaxation time scales or relaxation rates of all the vari- 
ables in the following sense. Let us characterize the long- 
time decay of <f>y(t) in the liquid by some time r„ . For 
example, as is occasionally done in analyzing molecular- 
dynamics simulation data, it may be defined as the center 
of the a process, i.e. as the time at which the correlator 
has decayed to 50% of its plateau value: 



fy(r Q 



fy/2. 



(19) 



finds 

<MCV) 



Y 



This time diverges upon approaching the glass transi- 
tion: in the leading asymptotic limit for a — ► —0, one 
Cyt'g.. Here, the constant Cy is defined by 
fy/2. All times or rates are proportional to 
each other, and follow a power law specified by the ex- 
ponent 7 defined in Eq. (|lq): 



l/rJ = IY 



(20) 



The constants of proportionality Cy or Ty depend on the 
variable Y and on the precise convention for the definition 



of , such as Eq. ( |19| ) . All the time scales or rates are 
coupled in the sense that the ratio Cy / Cy or Ty /Ty for 
two different variables Y and Y' becomes independent of 
the control parameters in the limit a — > —0. 

Figure [j] shows that for large elongations of the 
molecules, say £ > 0.6, and for very small elongations, 
say £ < 0.2, the exponent parameter A is close to the 
value 0.736 ± 0.003 for the HSS (C = 0). For C ~ Co 
the two contributions to the structure-factor peak for 
angular-momentum indices £ = and 1=2 are of equal 
importance jlj. For such elongations, A increases con- 
siderably. Unfortunately, Eq. (|8j) is so involved that we 
cannot trace back the increase of A to the underlying 
variations of the structure factor. 

The accurate determination of the time scale to is cum- 
bersome. It is given by the constant plateau value of 
the function = {[^ c {t) - /^ c ]Af } 1/o t within the 
time interval where the leading-order Eq. (11) is valid. 
The leading corrections to this law vary proportional to 
(to/t) 2a and cause deviations from the plateau at small 
times. The unavoidable errors in the determination of 
the critical value </? c (C) cause a ^ in Eq. (9b). Thus, 
^(t) increases proportional to t or —t for t — ► oo if a > 
or a < 0, respectively. For example, to determine to for 
C = 1.0 with an error estimated to be smaller than 3%, 
we have fixed y c (C) with 9 relevant digits. The plateau 
regime of $?(t) is largest for q = 9.8, where fg C has an 
intermediate value; it extends from t ~ 10 4 to t ~ 10 9 . 
We have checked that results for q = 3.0, 7.4, 13.0, and 
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FIG. 1: Exponent parameter A as function of the elongation 
£. The arrow indicates the critical elongation f e = 0.345 
separating the two glass phases. The dashed horizontal line 
marks the value A = 0.736 for the hard-sphere system. 



TABLE I: Parameters characterizing the MCT-liquid-glass- 
transition dynamics for systems with elongations £ = 0.0, 0.4, 
and 1.0. 



c 




C 


to 


A 


a 


b 


7 


B 


i- 


0.0 


0.530 


1.54 


0.0220 


0.736 


0.311 


0.582 


2.46 


0.838 


0.703 


0.4 


0.675 


1.81 


0.0123 


0.885 


0.222 


0.330 


3.77 


2.54 


0.110 


1.0 


0.565 


1.90 


0.0139 


0.739 


0.310 


0.576 


2.48 


0.857 


0.687 



16.2 lead to the same to within the specified error. A 
similar statement holds for the determination of to for 
other elongations. 

Table I compiles the parameters characterizing the uni- 
versal formulas for the three elongations to be considered. 



III. RESULTS FOR THE STRUCTURAL 
RELAXATION 

A. Density correlators 

Figures @ and || demonstrate the coherent density 



cor- 



(t) and the tagged-molecule correlation func- 



relators 

tions (f>q tS {t) and (j>g S (t) near the liquid-glass transition 



for the elongations C = 0.4 and 1.0, respectively. The 
results are for two representative wave numbers q: the 
wave number q « 7.0 is close to the first peak, and 
q 9.8 is near the first minimum of Sg (cf. Fig. 2 of 
Ref. ^) . The oscillatory transient dynamics occurs within 
the short-time window t < 1. The control-parameter sen- 
sitive glassy dynamics occurs for longer times for packing 
fractions <p near tp c . At the transition point <p — tp c , the 
correlators decrease in a stretched manner towards the 
plateau values (f^ , fgg or f^f) as shown by the dot- 
ted lines. Increasing tp above ip c , the long-time limits 
increase, as shown for the coherent correlators (j>^ (t) for 
q = 9.8. Decreasing ip below (f c , the correlators cross 
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^ = 0.4 




FIG. 2: Coherent density correlators cj>^ (t) [solid lines in 
(a) and (b)], tagged-molecule density correlators 4>^ a (t) [solid 
lines in (c) and (d)], and tagged-molecule charge-density cor- 
relators 4>q, a (t) [dashed lines in (c) and (d)] for the elon- 
gation £ = 0.4 as function of log 10 1 for two intermediate 
wave numbers q — 7.0 and 9.8. The decay curves at the 
critical packing fraction ip c are shown as dotted lines and 
marked by N and Z for the total density and charge den- 
sity correlators, respectively. Glass curves (e > 0) are shown 
only for cf>^ (t) with q = 9.8 in order to avoid overcrowding 
of the figure. The packing fractions are parameterized by 
e — (ip — tp c )/tp c — ±10~ x , and x = 1, 2, 3, 4 are chosen. 
The dashed-dotted line in (b) is a Kohlrausch-law fit for the 
x = 4 Q-process: f* c exp[-(t/r)' 3 ] with = 0.40. The filled 
circles and squares mark the characteristic times t a and t'^, 
respectively, defined in Eqs. (0) and (|l|) for x = 1, 2, 3 and 
4. Here and in the following figures, the diameter of the con- 
stituent atoms is chosen as the unit of length, d = 1, and the 
unit of time is chosen so that vt = 1. 



their plateaus at some time Tp, and then decay towards 
zero. For small | ip — (p c |, t@ is given by Eq. ([T^) for 
all the correlators. The decay from the plateau to zero 
is called the a-process. It is characterized, e.g., by the 
cv-time scale r a defined as in Eq. ( jl9| ) for each correla- 
tor. Thus, upon decreasing tp c — tp towards zero, the 
time scales and r Q increase towards infinity propor- 
tional to t a and t' a , respectively. The two-step-relaxation 
scenario emerges, because the ratio of the scales T a / t@ in- 
creases as well. Figures || and || exemplify the standard 




FIG. 3: Results as in Fig. [], but for the elongation £ = 1.0 
for wave numbers q — 7.4 and 9.8, and a Kohlrausch exponent 
(3 = 0.68. 



MCT-bifurcation scenario. For small | ip — ip c \, the re- 
sults can be described in terms of scaling laws mentioned 
in Sec. 



[IB. This was demonstrated comprehensively in 
Refs. |2j and || for the HSS, and the discussion shall not 
be repeated here. 

For £ = 1.0 and q > 5, the plateaus for the tagged 
molecule's total-density and charge-density fluctuations 
are very close to each other: f£? « f^° s (cf. Figs. 12 and 
13ofRef.0). Fi gures H(c) and 0(d) demonstrate that also 
the dynamics is nearly the same, 4>^ s (t) « (j)g S (t). This 
means that for q( > 5 and for strong steric hindrance, 
the cross correlations F^{t) are very small. The reason 
is that the intramolecular correlation factors jg(g£/2) are 
small, and thus interference effects between the density 
fluctuations of the two interaction sites are suppressed. 
Coherence effects can be expected only for small wave 
numbers. For this case, the functions can be understood 
in terms of their small-q asymptotes. The latter are de- 
termined by the dipole correlator and the mean-squared 
displacements and their results shall be discussed in 
the following two subsections. 

Figures ||(c) and ||(d) deal with tagged molecule's cor- 
relators of weak steric hindrance. In this case, the charge- 
density fluctuations behave quite differently from the to- 
tal density fluctuations. The most important origin of 
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this difference is the reduction of the mode-coupling ver- 
tex in Eq. (^) for x = Z relative to the one for x = N. 
For small elongations, the effective potentials from the 
surroundings for the tagged molecule's reorientation are 
small. Therefore, decreases strongly relative to /^f 
for C decreasing towards ( c (c/. Fig. 13 of Ref. |]). For 
C < (c, the charge-density fluctuations relax to zero as in 
a normal liquid. This implies, as precursor phenomenon, 
that the cv-time scale of the charge-density fluctua- 
tions shortens relative to the scale t* n for the total den- 
sity fluctuations. Thus, the differences between 4>q S (t) 
and <j>q s {t) for small £ shown in Figs. |^(c) and ||(d) are 
due to disturbances of the standard transition scenario 
by the nearby type- A transition from a normal glass to 
a plastic glass. 

The stretching of the relaxation process is much more 
pronounced for the C = 0-4 system than for the £ = 1.0 
system. The wave vectors q = 7.0 and 7.4 refer to the 
structure-factor peak position for £ = 0.4 and 1.0, re- 
spectively, and the corresponding plateau values are 
almost the same. Figure ||(a) demonstrates that the 
C = 0.4 correlator for e = — 10~ 4 requires a time in- 
crease by 2.3 decades for the decay from 90% to 10% of 
the plateau value fq C - The corresponding decay inter- 
val for the £ = 1.0 correlator is 1.7 decades, as shown 
in Fig. ||(a). Thus, the specified time interval for the 
a-process is 4 times larger for the small than for the 
large elongation. Often, stretching is quantified by the 
exponent (3y of the Kohlrausch-law fit for the a-process: 
(j)y{t > Tp) tx exp[— (t/r^)^}. Such fits are shown by the 
dashed-dotted lines for 4>q (t) for q — 9.8 and e = — 10~ 4 
in Figs. |(b) and |(b). For C = 0.4 and C = 1.0, the fit 
exponents (5 are 0.40 and 0.68, respectively, quantifying 
the larger stretching for the smaller elongation. The tran- 
sient dynamics for intermediate and large q is not very 
sensitive to changes of £. Both Figs. |^(b) and ||(b) for 
q — 9.8 demonstrate an initial decay of the correlators by 
about 20% if the time increases up to t w 0.1. This simi- 
larity is also reflected by the similarity of the microscopic 
time scale to for the two elongations listed in Table |. 
However, for e = — 10~ 3 or — 10~ 4 , the correlators <fiq s (t) 
require about 100 times longer times t for C, = 0.4 than 
for ( — 1.0 to decay to zero. All the enhanced stretching 
features reflect the fact that the anomalous exponents 
a and b are smaller for £ = 0.4 than for £ = 1-0 (cf. 
Table |). 

A comment concerning the accuracy of the numeri- 
cal solutions of the equations of motion might be in or- 
der. The primary work consists of calculating <fr^ (t) from 
Eqs. (||) and (0) on a grid of times. In the work reported 
here, the initial part of the grid consists of 100 values with 
the equal step size St — 10~ 5 . This interval is then ex- 
tended by successively doubling its length and the step 
size St. By inspection one checks | <f>^ (t) \ < 1 so that 

the Laplace transform (f>q (z) = i J °° dt exp(izf)^>f (t) 
exists as analytic function for Im z > 0. One checks 
for e < that <pq (t) decreases fast enough for large t so 

that ^(w) = lim^oO + %V) = tffH + *0f H 



-2 



-3 



1 1 

A 




i i i i i i i i 

S= 1.0, q = 7.4 


- / x=5 \ 


A 


log 10 X^"(to) _ 


7 


\ 


"V/ 
v./ 


■i i 


i i/ 


/"n 

i 



-12 -10 -8 -6 -4-2 2 



FIG. 4: Double logarithmic presentation of susceptibility 
spectra \q ( UJ ) as function of frequency uj for the density 
fluctuations discussed in Fig. ^(a) (see text). 

exists with a smooth reactive part <f>^ (uj) and a smooth 
non-negative spectrum cj>^ {uj). The Fourier integrals 
are evaluated with a simplified Filon procedure The 
solid lines in Fig. [| show results for the susceptibility 
spectrum \q ( w ) = w </>f ( w ) f° r C = 1-0, q = 7.4, 
and e = (if — <p c (.Q)/<Pc(Q = -10~ x with x = 4,5, 
and e = 0. The correlators are used to evaluate the 
polynomials mf (t) = J2k, P V q,kp(t>k (*)<^ (*) and from 
here one gets rriq (uj) = rriq (uj) + irriq (uj). Equa- 
tions (||) and (0) are solved if 4>q (w) agrees with (j>^ (uj) = 
-l/{u - (ttq) 2 /[ui + (Of ) 2 mf (w)]}. Because of causal- 
ity, it suffices to check 4>q (uj) = <pq (uj). The functions 

ujcjjq (uj) are shown as dotted lines in Fig. f|. The desired 
identity is verified by inspection within the accuracy bet- 
ter than that of the drawing. Corresponding statements 
hold for the other functions like <j>^' s z (i). The evalua- 
tion of (j)q (t) on the discrete set of times from Eqs. (^) 
and (0) is done by an extension of the method discussed 
in Ref. [l0| But it should be noticed that the described 
proof of the accuracy of the solution does not require an 
account of how the solution (j)q (t) has been obtained. 

B. Dipole correlators 

The dipole correlator of the tagged molecule is defined 

by 

C 1>s (t) = {e s (t)-e s (Q)). (21) 

Here, the unit vector e s denotes the tagged molecule's 
axis. Ci )S (t) is the reorientational-correlation function 
for angular- momentum index £ — 1. For a similar rea- 
soning as presented in the paragraph preceding Eq. (Q), 
it is identical to the coherent reorientational function. 
It can be obtained as the zero-wave-vector limit of the 
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FIG. 5: Dipole correlators Ci, s (i) for two elongations £ = 0.4 
(upper panel) and 1.0 (lower panel). The correlators at the 
critical packing fraction (p = <p c axe shown as dotted lines 
marked with c. The plateau ff a for each elongation is marked 
by a horizontal line. The distance parameters are chosen as 
e = (ifi — tp c )/(p c — — 10~ x with x — 2 (faster decay) and 
x = 3 (slower decay). The filled circles and squares mark the 
corresponding time scales t a and t' a , respectively, defined in 
Eqs. (0) and @. The open circles and squares on the curves 
mark the characteristic time scales r^'" and r^' s defined by 

CiA^'l = fts and Ci, a (ji' e ) = ftJ2, respectively. The 
vertical line for £ = 0.4 indicates the decay interval described 
by the asymptotic formulas for the /3 process (see text). 



charge correlator: </>f s (i) 



- C M (i)+0(g 2 ) §. C:.M) 
is evaluated most efficiently as follows. One carries out 
the q = limit in Eq. (|4|) for x = Z to get the exact 
Zwanzig-Mori equation 

d?c 1>s (t) + 2v%c hs (t) 

+ 2v 2 R f dt l m z s {t-t')dvC 1 , s {t l )=0, (22) 



to be solved with the initial condition C\. s (t) = 1 — 
(v R t) 2 + 0(\t\ 3 ). The kernel is the q = limit of the 
relaxation kernel m z s (t) from Eqs. (||): 



mf (t) = (C 2 /72^ 2 ) 



dfcfc 4 /< (cf) 2 ^ 



<t>ls(t)<f>k(t)- 

The integral is discretized to a sum over M terms as ex- 
plained in connection with Eq. (5b). Substituting the 



correlators </>% s (t) and <j>% (t) , the kernel mf (t) is deter- 
mined. The remaining linear integro-differential equa- 
tion (JH) is integrated to yield the desired result for 
Ci. s (t). Equation (|2"3|) yields directly the nonergodicity 



FIG. 6: Rescaled dipole correlators Ci, s (i) = [Ci, 3 (t) — 
fi,a]/hi,s (full lines) for two distance parameters e = (^ — 
tfic)/<Pc = — 1Q X with x = 2 (faster decay) and x = 3 (slower 
decay) . The dashed lines show the /3-correlators G(t) from 
Eq. (Q for the exponents A = 0.885 and 0.739 for the elon- 
gations £ = 0.4 and 1.0, respectively (cf. Table |). Here 
/i c , s = 0.69 and 0.97, and hi, s = 1.0 and 0.056 for ( = 0.4 and 
1.0, respectively. 



parameter of the kernel, /if = mf (t — ► oo), as integral 
over the products of f^ s fk- From /if, one derives the 
probability for the arrest of the dipole /i )S = C% >s (t — > 
) = /if /(l + /if )■ This number can also be obtained 



as /i . 



lim 9 ^o/ifs- The critical value /f s and the 



corresponding critical amplitude hi >s = linig->o ^ s were 
discussed in Fig. 13 of Ref. |l|. 

The dipole correlators Ci. s (t) for the elongations ( = 
0.4 and 1.0 are shown in Fig. [| for the critical point 
(p = ip c and for two liquid states. The time scales r^' s 
and t„' s characterizing the center of (5- and a-relaxation 
processes, denned by Ci )S (rg' s ) = /f s and Ci >s (r^ ,s ) = 
/f s /2, are marked by open circles and squares, respec- 
tively. The curves do not clearly exhibit the typical two- 
step-relaxation scenario. For £ = 1.0, the critical plateau 
is so large, f± a = 0.97, that only 3% of the decay are 
left for the transient motion and the critical relaxation. 
The results for £ = 0.4 are influenced by the nearby 
type- A transition. Let us consider the /3-relaxation pro- 
cess for the dipole dynamics. The factorization theorem, 
Eq. (|l2|), specializes to 



C ltS (t) = fl s + h ltS G(t). 



(24) 



This means that the rescaled correlators C\ >s {t) = 
[Ci )S (t) — /f iS ]/^i,s are given by the /3-correlator G(t) 
The latter obeys the scaling law, specified by Eq. (fi~3]) 
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by the q = limit of the mode-coupling functional s 
given in Eq. (6b) at the critical point. The a-time scale 



FIG. 7: Dipole correlators Ci, s (t) for the elongations £ = 0.4 
(upper panel) and 1.0 (lower panel) for various distance pa- 
rameters e — (ip — ip c )/ip c = — 10~ x with x specified in the fig- 
ure, presented as function of log 10 (t/T^' B ). The a-relaxation- 
time scale r*' 3 is defined by Ci, s (t„) = /f a /2. The horizontal 
lines indicate the plateaus /f s . 



and (14). For fixed rescaled time, t = t/t a , the cited 
formulas deal with the results correctly up to order yfa. 
The leading corrections are of order | a | , and they explain 
the range of validity of the leading results Q . Figure ^ 
demonstrates the test of the /3-scaling law. On a 10% ac- 
curacy level, the leading-order asymptotic law accounts 
for 12% of the decay of Cx tS (t) around the plateau for 
( = 0.4, while it describes only 1.4% of the decay for 
( = 1.0. The latter feature is due to the large critical 
nonergodicity parameter a , and hence, due to the small 
critical amplitude hi jS , for the strong steric hindrance. 
The decay interval described by the leading-order asymp- 
tote for £ = 0.4 is indicated in Fig. || by the vertical line. 
For e = —0.001, the corresponding dynamical window ex- 
tends from about t = 1.1 x 10 3 to about 2.4 x 10 4 , while 
it extends from about t = 20 to t — 70 for e = —0.01. 
This discussion requires a reservation. The corrections to 
the scaling results can lead to an offset of the plateau Q . 
One recognizes that incorporation of such offset would 
improve the agreement between numerical solution and 
its asymptotic description in the upper panel of Fig. |^. 
The mentioned windows are obtained only after a plateau 
offset had been eliminated. 

The a-relaxation scaling law for the dipole correlator 
reads according to Eq. (jig): 

Ci, s (t) = C hs (t), t = tli! a , k|«l, t a «t. (25) 

The e-independent master function Ci s is determined 



r^' s can be defined by Ci )S (t„ ,s ) — fi s /2, and shall be 
written as t^' s — ti^t' a . Here, the a-scale factor t\ tS is 
defined in terms of the master function as Ci, s (ii iS ) = 
fi s /2. The scaling law implies that a representation of 
C\ tS {t) as a function of the rescaled time t/r^ s should 
superimpose correlators for different distance parameters 
e on the common curve Ci ja (i/ii jS ). Asymptotic validity 
means that the log(i/T*' s ) interval, where the scaling law 
is obeyed, expands to an arbitrary size for e — > 0. 

The lower panel of Fig. demonstrates that the de- 
scribed scenario for the evolution of the a process is valid 
for the elongation £ = 1.0. On the other hand, the up- 
per panel shows that the dipole correlators for £ = 0.4 
exhibit the superposition principle only for | e | < 10~ 4 . 
In particular, the plateau region emerges only for these 
extremely small values of the distance parameter | e | . As 
discussed in Ref. this is because of the correction to 
the leading-order asymptotic law, which reads for not too 
large values of rescaled time t: 



C 1 , s (t)=C liS (t) + \a\B 1 h 1>s /? 



(26) 



with B x = - b)T(l + b) - A]. Therefore, the cor- 

rection is larger, the larger the product Bihi, s is. As 
discussed in connection with Fig. 13 of Ref. [j], the size of 
the critical amplitude /ii jS becomes very large near the 
critical elongation £ c . In addition, B\ increases with A. 
One gets B x = 1.56 and 0.45 for C = 0.4 and 1.0, re- 
spectively. Thus, the anomaly shown in the upper panel 
of Fig. is due to the large correction to the leading- 
order asymptotic law, caused by the precursor effect of 
the nearby type-A transition and by the large value for 
A. The descriptions of correlators by the first and second 
scaling laws overlap for t k, t a . Therefore, the anomaly 
for ( = 0.4 exhibited in Fig. |^ can also be explained as 
the large percentage of the decay of C\ jS (t) described by 
the /3-scaling law (cf. Fig. ||). 



C. Mean-squared displacements 

There are two mean-squared displacements (MSD) to 
be considered for the symmetric dumbbell. One refers 
to the position of the constituent atom of the tagged 
molecule r s and the other to the molecule's center r s c = 
(r s A + rf)/2: 



A A ,c(t) = ([r^(t)-f s ^(0)r)/6 



(27) 



Here, a factor 6 is introduced in the definition for later 
convenience. Since there is the relation || 

A A (t) = A c (t) + (C 2 /12) [1 - Ci,.(t)], (28) 

it is sufficient to calculate Ac(t). This function is de- 
termined by the small-g limit of the density correlator: 
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= 1 - 9 2 Ac(t) + 0(g 4 ) @. The small-g expan- 
sion of Eq. §Q) for x = N leads to the exact equation of 
motion 

d 2 t A c {t)-v 2 T + v 2 T / dt'm*(t-t')dt<A G (t') = 0, (29) 
Jo 

to be solved with the initial behavior Ac(t) = 
(l/2)(w T <) 2 + 0(|i| 3 )- The kernel is obtained as the q = 
limit of Eq. (|6a|) for x = N: 

poo 

m? (t) = (1/67T 2 ) / dfcfc 4 p^(cf) 2 W f < s (t)^fW. 
Jo 

After the discretization of the integral as explained above 
and substituting the results for the two density corre- 
lators, the kernel (t) is determined. Then, the lin- 
ear equation (|2^) is integrated to yield the MSD of the 
molecule's center. 

The long-time behavior of the MSD depends sensitively 
on control parameters near the transition singularity, and 
therefore it is well suited to study the glass-transition 
precursors. In liquid states, one obtains from Eq. j29| ) 
the long-time asymptote: lim^oo Ac (t)/t = D. Here 
D is the diffusion constant, and it is expressed as the 
inverse of the zero-frequency spectrum of the relaxation 
kernel: D = 1/ J Q °° dtm^ (t). In glass states, the tagged 
molecule's total density fluctuations arrest for long times: 
— * °°) = fqs > 0- The Lamb-Mossbauer factor 
fq 8 approaches unity for q tending to zero, and a local- 
ization length, rc , of the center characterizes the width 
of the /^,-versus-g curve: f^ a = 1 — (qrc) 2 + 0(q 4 ). 
One gets limt_ ) . 00 Ac(t) = r c . Using Eq. (p9|), one can 
express r 2 c as inverse of the long-time limit of the relax- 
ation kernel from Eq. (|3^): r 2 c — l/mf (t — > oo). 

The ideal liquid-glass transition implies a transition 
from a regime with molecule's diffusion for ip < tp c to 
one with molecule's localization for ip > cp c . The former 
is characterized by D > and 1/rc = 0, and the latter 
by D = and 1/rc > 0. The subtleties of the glass- 
transition dynamics occur outside the transient regime. 
Figure || exhibits the MSD of the molecule's center for 
various packing fractions near the transition. For very 
short times, say t < to, interaction effects are unimpor- 
tant and \im t ^o Ac (t) / 1 2 — reflects ballistic mo- 
tion. For times larger than i , the cage effect leads to 
a suppression of Ac (t) below the short-time asymptote. 
For long times in liquid states, the MSD approaches the 
diffusion asymptote, lim^oo Ac(t)/t = D, as shown by 
the dotted straight lines drawn for the curves with labels 
x = 1 and x = 4. Upon increasing ip towards ip c , the 
diffusivity decreases towards zero. 

The curves with x = 1 for e > in Fig. || deal with 
the glass state ip = l.lp c . For this density, there is no 
obvious glassy dynamics. Rather, Ac(t) has approached 
its long-time limit r c after the oscillations have disap- 
peared for t«l. Decreasing tp towards ip c , the softening 
of the glass manifests itself by an increase of the local- 
ization length rc- At the transition point ip = (p c , the 




FIG. 8: Double logarithmic presentation of the mean-squared 
displacement for the center of mass Ac (t) for the elongations 
£ = 0.4 (upper panel) and 1.0 (lower panel). The dotted lines 
with label c refer to the critical packing fraction ip c , and the 
solid lines to e = (<p — ip c )/iPc = ±W~ X for x = 1, 2, 3, 4. The 
straight dashed line with slope 2 in each panel exhibits the 
ballistic asymptote (vrt) 2 /2. The straight dotted lines with 
slope 1 denote the long-time asymptotes Dt of the two liquid 
curves for x = 1 and 4. The horizontal lines mark the square 
of the critical localization lengths (r c ) 2 . The filled circles and 
squares mark the characteristic times t a and t' a , respectively, 
defined in Eqs. @ and @ for x = 1, 2, 3, 4. 

critical value r c c = 0.0587 (0.0497) for C = 0.4 (1.0) is 
reached. Substituting Eq. jlQ ) for f** and the analogous 
formula for f^ s to evaluate r c from Eq. (^fj|), it follows 
that the glass instability at ip c causes a -^cr-anomaly for 
the localization length, 

r 2 c = (r c f-h c V<7/(1-A) + 0(<t), (31) 

where he — 9.66 x 10~ 3 (5.14 x 10~ 3 ) for the elongation 
£ = 0.4 (1.0). A corresponding formula holds for the 
square of the localization length ta for the constituent 
atom with (r c f replaced by (r c A ) 2 = (r c ) 2 + (C 2 /12)(l - 
fi a ) and he replaced by hA = ^c + (C 2 /12)/ii iS according 
to Eq. (p8|). In Fig. [)| the leading order results for r c and 
r\ are shown as dashed and dotted lines, respectively. 
For ( — 1.0, the data are described by the square-root 
law for e < 2 x 10~ 3 . Similar small intervals for the 
validity of the leading order descriptions have been found 
for the Debye- Waller factor of the HSS for small wave 
vector q Q. For £ = 0.4, the range of validity of the 
universal formula is reduced to the even smaller intervals 
e < 0.5 x 10~ 3 . 
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FIG. 9: Square of the localization lengths for the molecule's 
center r% (circles) and for the constituent atom r\ (squares) 
for e = (<p- ip c )/tp c = 10~ n/3 , n = 7,8, ■ ■ ■ , 12. The crosses 
mark the critical values {tq) 2 and {r A ) 2 , respectively. The 
dashed and dotted l ines are th e leading order asymptotic laws 
r\ = ( r x) 2 ~ hx\/o-/(l — A) for X — C and A, respectively. 

The glass curves for e = 0.01, shown in Fig. || with 
label x = 2, exhibit a decay between the end of the tran- 
sient oscillations and the arrest at r 2 c which is stretched 
over a time interval of about two orders of magnitude. 
A similar two-decade interval is needed for the liquid 
curves (e < 0) with label x = 2 to reach the critical 
value (r^) 2 . After crossing (r£.) 2 , two further decades 
of an upward bent log 10 Ac (i)-versus-log 10 t variation is 
exhibited before the final diffusion asymptote is reached. 
The stretching is more enhanced as the system is driven 
towards the transition point, \e\ — > 0. The indicated 
slow and stretched time variation outside the transient 
regime is the glassy dynamics exhibited by the MSD. 
Again, the stretching is more pronounced for £ = 0.4 
than for C = 1.0. 

The first scaling law deals with the dynamics where 
Ac(t) — (r c ) 2 is small. One gets in analogy to Eqs. fl^ ) 
and (|l|): 

Ac(t) = (r c c ) 2 -hcVW\9±(t/U), a^0, (32) 

for | a | « 1 and t 3> to- For glass states with a > 0, 
this formula describes the approach towards the arrest at 
r 2 c = Ac{t — > oo). For liquid states with a < 0, it de- 
scribes how the Ac(£)-versus-i curve crosses and leaves 
the plateau (r^) 2 . In particular, one gets von Schwci- 
dler's law for large t/t a in analogy to Eq. (|l5|): 

AcW = {r c c ) 2 + h c (t/t'J, a - -0, t a « t « t' a . 

(33) 

The increase of Ap(£) above the plateau towards the dif- 
fusion asymptote is the a process of the MSD. In analogy 




FIG. 10: Double logarithmic presentation of the mean- 
squared displacement for the constituent atom AA(t) and for 
the center Ac(t) for the elongations £ = 0.4 (upper panel) 
and £ = 1.0 (lower panel). The curve A^(t) is shifted up- 
wards by two decades to avoid overcrowding. The distance 
parameter is e — — 10 -3 , and the corresponding times t a and 
t' a are marked by filled circles and squares, respectively. The 
dashed lines are the first-scaling-law results, Eq. (|3^). The 
open diamonds mark the points where the dashed lines dif- 
fer from the solid ones by 5%. The straight dashed-dotted 
lines exhibit the diffusion asymptotes Dt, and the filled di- 
amonds mark the position where these differ from the solid 
lines by 5%. The dotted lines exhibit the second-scaling-law 
results, Eq. (||). The horizontal lines mark the square of the 
localization lengths {r c A ) 2 and (r^) 2 . 

to Eq. (|l8|), there holds the superposition principle 
A c (t) = A c (t), i = t/t' a , | a | < 1, t„ < t. (34) 

In Ref. [ll], it has been discussed in detail how these 
leading-order asymptotic results can account quantita- 
tively for the glassy dynamics of the MSD, albeit for the 
tagged spherical particle and the dumbbell molecule im- 
mersed in the hard-sphere system. Instead of repeating 
such an extensive analysis, we shall only show here some 
representative examples. 

The upper and lower panels of Fig. [Io[ respectively for 
the elongations £ = 0.4 and 1.0, exhibit such tests of the 
scaling-law descriptions for a liquid state which is suffi- 
ciently close to the transition point, e = — 10~ 3 . Also, 
the MSD of the constituent atom A^(i) is considered 
in these figures. Because of Eq. (p8[), the asymptotic 
laws ( p2| ) and (J33) for A^(i) hold with (r^) 2 replaced 
by (r C AT - (r c c7+ (C 2 /12)(l - /f,J and h c replaced 
by Ha — he + (C 2 /12)/ii. s . The range of validity of the 
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first-scaling-law description is indicated by the open dia- 
monds. One finds that such a range is nearly the same for 
both Ac(t) and Aj±(t), and this holds irrespective of the 
molecule's elongation £. The beginning of the a-process 
of the MSD, i.e., its initial increase above the plateau, 
is described by von Schweidler's law, Eq. (|3^). It is ex- 
hibited for t ^> t„ by the dashed lines. The a-process 
terminates in the diffusion law for long times, exhibited 
by the straight dashed-dotted lines. The beginning of the 
diffusion law is indicated by the filled diamonds for each 
MSD curve. The a-process follows well the second scaling 
law, Eq. (HJ), which is presented by the dotted lines. The 
descriptions by the two scaling laws overlap for t » t a . 
Together, the two scaling laws provide a description of 
the glassy dynamics of the MSD for t > 100 and t > 10 
for £ = 0.4 and £ = 1.0, respectively. Notice that there 
is a large interval of times outside the transient regime, 
say log 10 1 > —0.5, where the structural relaxation is not 
described by the first scaling law. This is a peculiarity of 
the MSD [jll| which is not found for the other functions 
discussed here. 

The crossover window from the end of the von- 
Schweidler-law regime (indicated by the right open di- 
amonds) to the beginning of the diffusion process (indi- 
cated by the filled diamonds) for the MSD of the center, 
A c (t), is about 3 (2) decades wide for C = 0.4 (1.0). On 
the other hand, the corresponding window for the MSD 
of the constituent atom, A^t), is larger than that of 
A c (t) by a factor of about 7 (13) for C = 0.4 (1.0). As 
discussed in Ref. |ll|, this enlarged crossover window for 
A^(t) is caused by the rotation-translation coupling, i.e., 
by the second term on the right-hand side of Eq. (28). 
This effect is smaller for £ — 0.4 because the relaxation 
of the second term, determined by that of Ci yS (t), is con- 
siderably enhanced due to the nearby type- A transition 
as discussed in connection with Fig. o. 



Figure 11 exhibits the critical localization lengths of 



the molecule's center r c c and of the constituent atom r A 
as function of the elongation £. There are three distinct 
regions: one for small elongations, say £ < 0.3, another 
for large elongations, say £ > 0.4, and the crossover re- 
gion between them. For £ < £ c , the localization length 
of the center r c c is nearly constant while that of the con- 
stituent r c A increases rapidly with £. Both r£. and r° A 
decrease rapidly within the crossover region, and they 
decrease only slightly as function of the elongation for 
£ > 0.4. The rapid increase of r c A for £ < £ c can 
be understood as follows. Since C\ yS (t — > oo) = for 
£ < Co one gets from the long-time limit of Eq. (28): 
r c A = ^J{r c c ) 2 + C 2 /12- The dotted line in Fig. [ll] shows 
this formula with r c c fixed to the value for £ = 0.0. It ex- 
plains the increase of r c A for £ increasing up to £ c . Thus, 
r c A increases because the reorientational dynamics of the 
molecule's axis is not arrested for £ < £ c . In this sense, 
the localization of molecules with small elongations is 
primarily caused by that of the molecule's center. This 
is consistent with the result that the glass transition for 
systems with small elongations is mainly driven by the 
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FIG. 11: Critical localization lengths for the constituent 
atom r c A (solid line) and the center r c c (dashed line) along the 
liquid-glass-transition line parameterized by the elongation 
£. The dotted line represents the function yj (r^) 2 + £ 2 /12 
discussed in the text. The inset exhibits 1 — /f, s as function 
of £, where /f s = Ci lS (t — > oo) is the long-time limit of 
the dipole correlator Ci jS (t). The arrow indicates the critical 
elongation £ c = 0.345. 



arrest of the center-of-mass density fluctuations |l| . The 
critical localization length r c c provides the upper limit 
for rc characterizing the size of the arrested structure, 
and its value 0.07 ~ 0.08 for £ < £ c is consistent with 
Lindcmann's melting criterion. The critical localization 
length of the constituent atom r c A for £ > 0.4 has similar 
values. This implies that the localization of molecules 
with large elongations is caused by that of the molecule's 
constituents rather than by that of the molecule's center: 
the localization of the center is subordinate to that of the 
constituents. As explained in Ref. [l| angular correlations 
become more relevant for the glass transition of systems 
with large elongations. Such angular correlations result 
in more efficient localization of the constituent atoms, 
and this is the reason for the relevance of r c A for the local- 
ization of molecules with large elongations. The center 
localization is reduced by about l/v^, as expected for 
independent motion of the two constituents. 



D. The a-relaxation scales 

The superposition principles for various correlators im- 
ply coupling of the a-relaxation time scales or relaxation 
rates as described in Sec. [IB. This scale coupling or a- 
scale universality is demonstrated in Fig. 



J for the rate 

of the coherent total density correlator <fi^ (t) for 
the wave number q « 7, the rate l/r^ ,s of the dipole cor- 
relator Ci jS (£), and the diffusion constant D. Here, the 
a-relaxation times for 4>q{t) and Ci jS (t) are determined 
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FIG. 12: Double logarithmic presentation of the a-relaxation 
rates l/r a (squares) for the coherent total density correlator 
4>q {t) at the structure-factor-peak position, 1/t*' s (circles) 
for the dipole correlator Ci, s (i), and the diffusion coefficients 
D (diamonds) for the elongations £ = 0.4 (upper panel) and 
1.0 (lower panel) as function of the reduced packing fraction 
[ip c — <p)/<p c = 10~"/ 3 . The a-scaling times t& and r„' s 
are defined as in Eq. (ha). The solid lines are the power- 
law asymptotes Ta\ e | 7 (see text) with 7 = 3.77 and 2.48 for 
f = 0.4 and 1.0, respectively. The prefactors Ta have been 
chosen so that the solid lines go through the data points for 
n = 12. 




_4 I 1 1 1 1 1 1 1 

-2 -1.5 -1 -0.5 

log 10 ((P c -<P)/<Pc 

FIG. 13: Results as in Fig. but only for a restricted region 
of the reduced packing fraction, n < 6, for the elongation 
£ = 0.4. The solid lines here are the power-law fits T c jf e 7 
with the effective power-law exponents 7 eff = 2.90, 2.50, and 
2.05 for 1/t„ , D, and 1/tq' s , respectively. The prefactors 
T c a have been chosen so that the solid lines go through the 
data points for n = 5. 



by the convention given in Eq. (|19|). 

Let us first consider the results for £ = 1.0. Although 
the asymptotic behavior is the same for all the quan- 
tities, 1/t^' s and D start to deviate visibly from their 
asymptotic results for n = 4, while starts to devi- 

ate only for n = 2. The mentioned a-scale universality 
holds in the leading asymptotic limit for a — > — 0, and 
the corrections to the asymptotic predictions are differ- 
ent for different quantities. Thus, the found feature in 
the results for £ = 1.0 underlines the nonuniversality of 
the deviations. Let us add that the range of validity of 
the asymptotic-law description is quite similar to that of 
the hard-sphere system 0, |). The results for £ = 1.0 
follow the pattern that has been analyzed theoretically 
so far. The results for Q — 0.4 exhibit much more pro- 
nounced deviations from the asymptotic-law predictions: 
1/t„ and D start to deviate visibly for n = 7, and the 



deviation of 1/t^ ,s starts even for n — 9. These devi- 
ations are due to the corrections to the cv-scaling law. 
Their magnitudes are proportional to the product of the 
critical amplitudes for each quantity and the coefficient 
B\ as noted in Eq. (|2"6|). The critical amplitude for 
C = 0.4 at the structure-factor-peak position is slightly 
larger than the corresponding result for ( = 1.0 (cf. Fig. 9 
of Ref. |l|). In addition, B\ is 3.5 times larger for £ = 0.4 
than for 1.0. This is the reason for the larger deviation 
of 1/t„ found for ( = 0.4 than the one for £ = 1.0, and 
similarly for D. The even more pronounced deviation for 
l/ri' s is due to the large critical amplitude hi iS caused 
by the precursor effect of the nearby type-^4 transition, 
as discussed in connection with Fig. 13 of Ref. |l|. An- 
other possible source for the large deviations for ( = 0.4 
might be due to higher-order glass transition singulari- 
ties, which can lead to the violation of the second scaling 
law, and thus the a-scale universality. The signature 
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FIG. 14: Products Dr" (squares) and -Dt,J ,s (circles) for the 
a-relaxation scales of the £ = 0.4 liquid shown in Figs, [li] and 
|l|. The full (dashed) line is a power-law fit T'| e | 7 ' with V 
chosen so that the line goes through the data point for n — 5 
and i = -0.40 (0.45). 



of such singularities is the approach of A towards unity, 
which implies a divergency of B\. 

A remark shall be added concerning the determination 
of the exponent 7 entering the power-law behavior for 
the a-relaxation time scale or relaxation rate as specified 
by Eq. (|20|). This result is based on the validity of the 
scaling law. Therefore, one cannot appeal to MCT if one 
fits power laws for a-relaxation rates for cases where the 
scaling law is violated so strongly as shown in the upper 
panel of Fig. [?] for x < 2, i.e. for n < 6. Figure p"3| 
demonstrates that, for £ = 0.4, the a-relaxation rates 
for I e I > 10 -2 can be fitted well by power laws 1/r oc 

off 

I e I 7 for an about 1.3-decade variation of the distance 
parameter |e|. The used effective exponent 7 cff = 2.90 
(2.50; 2.05) describes the variation of 1/r^ (D; I/t*' 8 ) 
over 3.5 (3.3; 2.5) orders of magnitude. The variable- 
dependent effective exponent 7 cff again underlines the 
nonuniversality of the deviations from the scaling law. 
The quantities 7 cff are mere fit parameters without well- 
defined meaning for the discussion of our theory. They 
depend on the interval of | e | chosen for the fit. 

Another manner for testing the coupling between the 
scales of two variables, say Y and Y', is provided by a plot 
of the ratio /t£ versus the control parameter. Such a 
plot does not require the knowledge of co c , nor is it biased 
by the choice of some fit interval. This is demonstrated 
in Fig. |lj for the ratios formed with or t^' s and the 
scale for the diffusivity oc 1/D. Instead of weakly 
co-dependent ratios expected for the asymptotic law, the 
figure shows variations by more than a factor 3. The 
shown plot suggests power-law fits t% /t% oc (co c — ip)" 1 as 
shown by the full line with 7' = —0.40 and the dashed one 
with 7' = 0.45. The fit exponents 7' are the differences 
of the effective exponents 7 eff referring to t£ and t£ . 



IV. CONCLUSIONS 

The recently developed mode-coupling theory (MCT) 
for molecular systems has been applied to calculate 
the standard correlation functions which demonstrate 
the evolution of glassy dynamics of a symmetric-hard- 
dumbbell system (HDS). The equilibrium structure of 
this system is determined solely by excluded-volume ef- 
fects, i.e., the sluggish dynamics and the structural arrest 
are due to steric hindrance for translational and reorien- 
tational motions. The phase diagram demonstrates that 
there are two scenarios for the liquid-glass-transition dy- 
namics: one deals with strong and the other with weak 
steric hindrance for reorientational motion |l| . 

For the strong-steric-hindrance scenario, as it is ob- 
tained for dumbbells of elongation £ exceeding, say, 0.6, 
the constituent atoms are localized in cages of a similar 
size as found for the simple hard-sphere system (HSS) 
and for the motion of a single dumbbell in the HSS, 
Fig. [ll]. The parameter A, which determines the anoma- 
lous exponents of the decay laws and the critical time 
scales, is close to that of the HSS, Fig. [l]. The range of 
validity of the scaling laws for the (3- and a-processes, 
exemplified in the lower panels of Figs. |^ and [^.respec- 
tively, is similar to what was found for the HSS |H, jl and 
for the motion of a hard dumbbell in the HSS || |TT| . 
The rotation-translation coupling implies for the mean- 
squared displacement of the constituent atom a crossover 
interval between the end of the von Schweidler process 
and the beginning of the diffusion process which is about 
an order of magnitude larger than for the motion of a 
sphere in the HSS. This effect is even more pronounced 
for the dumbbell liquid, the lower panel of Fig. [l(], than 
for a dumbbell moving in the HSS tyl. Testing this pre- 
diction by a molecular-dynamics simulation would pro- 
vide a valuable information on the relevance of our the- 
ory. 

The results of this paper together with the preced- 
ing findings on the a-peaks for reorientational motion 
for angular-momentum index t = 1 and I = 2 as well as 
for the elastic modulus Q lead to the conclusion that the 
strong-steric-hindrance scenario explains the qualitative 
features of the structural relaxation in glass-forming van- 
der-Waals systems like orthoterphenyl, Salol, or propy- 
lene carbonate. This holds with two reservations. First, 
the calculated wave-vector dependence of the Debye- 
Waller factor is stronger than that measured by 
neutron-scattering spectroscopy for orthoterphenyl [ jT2| . 
It remains to be shown that application of the theory to 
molecules more complicated than dumbbells will lead to 
a smearing out of the strong q- variations characteristic of 
hard-sphere-like systems. Second, molecules which allow 
for a measurement of the dipole-correlator by dielectric 
loss spectroscopy carry an electric dipole moment. This 
leads to long-ranged interactions between the molecules. 
It remains to be studied how the incorporation of these 
interactions change the results based solely on hard-core 
interactions. 
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The dipole-susceptibility spectra for a dumbbell with 
large £ moving in a HSS || obey the scaling proposed by 
Dixon et al. jlj|. We found this is the case also for the 
spectra of the HDS for elongations in the range 0.6 < £ < 
0.8. For £ = 1.0, on the other hand, the spectra for the 
high-frequency wing are below the master curve found in 
Rcf. 13, and for ( — 0.6 they are slightly above. Thus, 
our theory is not consistent with the conjecture that the 
quoted scaling law holds universally. 

Decreasing the elongation, the steric-hindrance effects 
for reorientations weaken. A new scenario emerges for 
sufficiently small £ for two reasons. First, there is a criti- 
cal elongation £ c separating a normal glass from a plastic 
one. For ca = (C — Cc)/Cc approaching zero, the criti- 
cal amplitudes h^ s , in particular /ii jS = lim g __»o hg [ s , di- 
verge. This transition is due to the top-down symmetry 
of the molecules. But the asymptotic formulas for small 
€a remain valid also for nearly symmetric molecules; one 
merely has to replace | | by y/ e A + S 2 , where S is a 
number quantifying the symmetry breaking jl4|, [l^] . The 
precursor phenomena of the approach of £ towards Q c 
have been discussed before for the motion of a single 
dumbbell in a HSS [§). The most obvious one is the 
speeding up of the a-process for reorientations relative to 
that for density fluctuations, as demonstrated in Fig. [|. 
The second reason is the strong increase of the exponent 
parameter A for £ ~ ( Cl Fig. [|. For A = 1, there oc- 
curs a higher-order glass transition singularity, where the 
asymptotic expansions of the MCT solutions are utterly 
different from those cited in Sec. [IB. The approach of A 



towards unity leads to a shrinking of the range of validity 
of the discussed universal formulas. The corrections di- 
verge, as was discussed in connection with the coefficient 
B\ in Eq. (^(|). This equation shows that the increase 
of the product B\hi^ s leads to an increasing violation of 
the superposition principle for the a-process, as shown in 
the upper panel of Fig. [j| As a result, the a-relaxation 
time approaches the universal law, Eq. (^0|), only for very 
small distance parameters e = (ip — (/j c (£))/<^ c (£). For 
C = 1-0, the power-law regime is reached for | e | = 10 _1 , 
while | e | has to be smaller than 10~ 2 for the approach 
to the asymptotic law for ( = 0.4, Fig. [l^. A decrease of 
| e | by an order of magnitude is equivalent to an increase 
of the relaxation time by more than a factor of 1000. 

For £ = 0.4, the variation of the a-relaxation-time 
scales over more than three orders of magnitude can 
be described well by power laws, Fig. [l3|. But, the 
above-explained corrections to the leading-order asymp- 
totic laws imply that the fitted effective exponents 7° ff 
are considerably smaller than the exponent 7 = 3.77 
specifying the correct asymptotic behavior. Moreover, 
the corrections depend on the variable considered and 
therefore the found effective exponents are pairwise dif- 
ferent. This means that the a-relaxation scales are not 
coupled, rather the ratio of two scales varies according to 
a power law as shown for two cases in Fig. [l4|. Molecular- 
dynamics-simulation studies for a binary Lennard-Jones 
system which have been started by Kob and Ander- 



sen [[16| have been used for detailed tests of MCT, as can 
be inferred from Ref. [l?] and the papers cited there. The 
range of time variations for the a-process in these studies 
is about the one considered in Fig. As a remarkable 
deviation from the a-scale universality, a deviation of the 
fit exponent for the diffusivity from the exponent derived 
from density- relaxation curves has been reported p6[ . 
Indeed, the Dr Q -diagram for the simulation data is in 
semi-quantitative agreement with the results shown in 
Fig. |lj |lq| . Therefore, it is tempting to conjecture that 
the simulation data can be explained as done above for 
the results in Figs. |l3| and [l4|. 

We are not aware of experiments which exhibit the 
weak steric hindrance scenario. Molecular-dynamics sim- 
ulations for a liquid of symmetric Lennard-Jones dumb- 
bells with elongation £ = 0.33 lead to the suggestion |ll| 
that an inclusion of the rotational degrees of freedom is 
decisive for an understanding of the exponent parame- 
ter. This conclusion is corroborated by Fig. [I]. Detailed 
simulation studies of the glassy dynamics of a liquid of 
slightly asymmetric Lennard-Jones dumbbells for £ = 0.5 
have been reported by Kammerer et al. j2(], |2l], |2^] . The 
correlation functions dealing with the translational de- 
grees of freedom for large and intermediate wave vectors 
and also the ones for reorientational fluctuations for large 
angular-momentum indices could be interpreted with the 
universal MCT formulas. However, the dipole correlator 
did not show a two-step-relaxation scenario and it exhib- 
ited strong violations of the superposition principle quite 
similar to what is shown in the upper panel of Fig. [7]. 
The a-relaxation time scales for density fluctuations of 
intermediate wave vectors, for the diffusivity, and for the 
dipole correlator could be fitted well by power laws with 
the exponents 2.56, 2.20, and 1.66, respectively. The 
differences in these exponents are quite similar to what 
is demonstrated in Fig. [l3] for the corresponding quan- 
tities. This explains why the Dr a -versus-temperature 
diagram for the simulation results (2^] shows violations 
of the scale coupling in semi-quantitative agreement with 
the ones shown in Fig. [l4|. It seems that the simulation 
results for Ci jS (i) also fit nicely into the framework of the 
ideal MCT for molecular liquids. 

Let us notice that the density correlators of the glass 
states for wave vector q = 9.8 exhibit oscillations for 
times around 0.1, Figs. || and |[ These are the ana- 
logues of the oscillations analyzed previously for the HSS 
in connection with a discussion of the so-called boson- 
peak phenomenon and high-frequency sound |24|]. It 
should be mentioned that the dynamics of a dipolar- hard- 
sphere system was analyzed recently within the mode- 
coupling theory |25|] describing the structure by tensor- 
density fluctuations. Some comments on the general rela- 
tion between this theory and the one used in the present 
paper can be found in Rcfs. [I] and |2^. For the dipolar- 
hard-sphere system, the oscillations have been analyzed 
in detail. They reflect subtle couplings between transla- 
tional and rotational degrees of freedom as enforced by 
the long-ranged Coulomb interactions. For the HDS such 
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couplings are present as well, albeit caused by the short- and R. Schilling for helpful critique of the manuscript, 
ranged steric hindrance effects; but an analysis of these 
oscillations remains to be done. 
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